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We investigate in this paper properties of the spin resistivity in the magnetic semiconducting 
MnTe of NiAs structure. MnTe is a crossroad semiconductor with a large band gap. It is an an- 
tiferromagnet with the Neel temperature around 310K. Due to this high Neel temperature, there 
are many applications using its magnetic properties. The method we use here is the Monte Carlo 
simulation in which we take into account the interaction between itinerant spins and lattice Mn 
spins. Our results show a very good agreement with experiments on the shape of the spin resistivity 
near the Neel temperature. 
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I. INTRODUCTION 

Spin resistivity in materials has been a subject of in- 
tensive studies both experimentally and theoretically for 
more than five decades. Experiments have been perfor- 
med to determine the spin resistivity p in many magne- 
tic materials from metals to semiconductors. The rapid 
development of the field is due mainly to many appli- 
cations in particular in spintronics. One interesting as- 
pect of magnetic materials is the existence of a magne- 
tic phase transition from a magnetically ordered phase 
to the paramagnetic (disordered) state. Depending on 
the material, p can show a sharp peak at the magne- 
tic transition temperature Tc/ or just only a change 
of its slope, or an inflexion point. The latter case gives 
rise to a peak of the differential resistivity dp/dT^ 
Very recent experiments such as those performed on 
ferromagnetic SrRuOa thin films^, Ru-doped indu- 
ced ferromagnetic Lao.4Cao.6Mn03^, antiferromagnetic 
e-(Mni_a^Fea;)3.25Ge^, semiconducting Pro.rCao.sMnOa 
thin films^, superconducting BaFe2As2 single crystals^, 
and Lai-a^Sr^^MnOa^ compounds show different forms of 
anomaly of the magnetic resistivity at the magnetic phase 
transition temperature. 

The properties of the total resistivity stem from dif- 
ferent kinds of diffusion processes : the scattering of the 
itinerant spins by phonons, by lattice magnons, by impu- 
rities and defects etc. Each contribution has in general a 
different temperature dependence. Let us summarize the 
most important contributions at low temperature (T) in 
the following expression 



p{T) = po + AT^ + + Cln ^ 



(1) 



where A, B and C are constants. The first term is T- 
independent, the second term proportional to repre- 
sents the scattering of itinerant spins at low T by lattice 
spin- waves. Note that the resistivity caused by a Fermi 
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liquid is also proportional to T^. The term corres- 
ponds to low-T resistivity in metals. This is due to the 
scattering of itinerant electrons by phonons. Note that at 
high T, metals show a linear-T dependence. The In term 
is the resistivity due to the quantum Kondo effect at very 
low T. For the magnetic contribution to the total resis- 
tivity, the term has been obtained from the magnon 
scattering by Kasuya^^. However, at high T in particular 
in the region of the phase transition, much less has been 
known. The general idea that the magnetic resistivity is 
a function of the spin-spin correlation was introduced by 
de Gennes and Friedel.^^ According to this idea, the ma- 
gnetic resistivity should behave as the magnetic suscepti- 
bility, thus it should diverge at T^. Fisher and Langer— , 
and Kataoka^^ have suggested that the range of spin-spin 
correlation changes the shape of p near the phase transi- 
tion. The resistivity due to magnetic impurities has been 
calculated by Zarand et al.^^ as a function of the An- 
derson's localization length. This parameter expresses in 
fact a kind the correlation sphere induced around each 
impurity. Their result that the resistivity peak depends 
on this parameter is in agreement with the spin-spin cor- 
relation idea. In our previous works ^^^^ we have studied 
the spin current in ferromagnetic thin films by Monte 
Carlo (MC) simulations. The behavior of the spin resisti- 
vity as a function of T has been shown to be in agreement 
with main experimental features and theoretical investi- 
gations mentioned above. We have introduced in these 
works the picture of scattering of itinerant spins by ma- 
gnetic defect clusters which are known to be formed in 
the transition temperature region. The size of each defect 
cluster expresses a kind of correlation between spins. 

In antiferromagnets much less is known because there 
have been very few theoretical investigations which have 
been carried out. Haa^ has shown that while in fer- 
romagnets the resistivity p shows a sharp peak at the 
magnetic transition of the lattice spins, in antiferroma- 
gnets there is no such a peak. The alternate change of 
sign of the spin-spin correlation with distance may have 
something to do with the absence of a sharp peak. 

In this paper, we are interested in the antiferromagne- 
tic MnTe, a well-studied semiconductor with numerous 
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applications due to its high Neel temperature. The pure 
MnTe has either the zinc-blend structure^^ or the hexago- 
nal NiAs one^^. We confine ourselves in the case of hexa- 
gonal structure. For this case, the Neel temperature is 
Tat = 310 K^^. Hexagonal MnTe is a crossroad semicon- 
ductor with a big gap (1.27 eV) and a room-temperature 
carrier concentration of n = 4.3x lO^^cm^^^^s^^^ Without 
doping, MnTe is non degenerate. In doped cases^^^^, 
band tails created by doped impurities can cover more or 
less the gap. But these systems, which are disordered by 
doping, are not a purpose of our present study. So, in the 
following we study only the pure MnTe. The behavior of 
the spin resistivity p in MnTe as a function of T has been 
experimentally shown^~— . In our previous paper, we 
have theoretically studied with the Boltzmann's equation 
using numerical data for cluster distribution obtained by 
MC simulations. We could only compare our result with 
old experimental data available at that time below the 
transition temperature.^^ After the publication of our 
paper, new experimental results become available for the 
whole range of temperature for pure MnTe performed by 
He et al^^ This motivates the present work. 

In this paper, we use the same model as in our pre- 
vious work^^ for MnTe. But unlike that previous paper 
which used an approximate Boltzmann's equation, we 
shall use here direct MC simulations without recourse 
to approximations. In addition, we take into account the 
lattice magnetic relaxation time in the calculation of the 
spin resistivity. The tuning of the relaxation time allows 
an excellent agreement with experimental resistivity of 
MnTe as will be seen below. 

In section [Hi we show our model and describe our MC 
method. Results are shown and compared with experi- 
mental resistivity of MnTe in section IIIII Concluding re- 
marks are given in section HVl 




Figure 1: Structure of MnTe of NiAs type is shown. Anti- 
parallel spins are shown by black and white circles. Nearest- 
neighbor (NN) interaction is marked by Ji, next NN interac- 
tion by J2, and third NN one by J3. 



II. MODEL 

The hexagonal NiAs-type structure of MnTe is shown 
in Fig. 1.^^ It is composed of ferromagnetic xy hexagonal 



planes antiferromagnetically stacked in the c direction. 
The nearest-neighbor (NN) distance in the c direction 
is c/2 3.36 A shorter than the in-plane NN distance 
which is a = 4.158A. Neutron scattering experiments 
show that the main exchange interactions between Mn 
spins in MnTe are i) interaction between NN along the 
c axis with the value Ji/ks = —21.5 ± 0.3 K, ii) ferro- 
magnetic exchange 321^3 ^ 0.67 ±0.05 between in-plane 
neighboring Mn (they are next NN by distance), iii) third 
NN antiferromagnetic interaction J^/ks — —2.87 ±0.04 
K. 

The interactions Ji, J2 and J3 are indicated in Fig. [TJ 
In addition, the spins are lying in the xy planes perpen- 
dicular to the c direction with a small in-plane easy-axis 
anisotropy D?^ We note that the values of the exchange 
integrals given above have been deduced from experimen- 
tal data by fitting with a formula obtained from a free 
spin- wave theory Other fittings with mean-field theo- 
ries give slightly different values. Therefore, care should 
be taken while using values deduced from experimental 
data, keeping in mind that they depend on models used 
and approximations involved in the fitting. 

The lattice Hamiltonian is given by 

1-L = — Ji ^ S^.Sj — J2 ^ Si.Sm — J?, ^ S^.S/c 

(i,m) (i,/c) 

-DY,{S^? (2) 

i 

where is the Heisenberg spin at the lattice site ^, J],^^ ^-^ 
is made over the NN spin pairs and Sj with interaction 
Ji, while ^(^^ ^) and ^^^^ are made over the NNN and 
third NN neighbor pairs with interactions J2 and J3, res- 
pectively. I) > is an anisotropy constant which favors 
the in-plane x easy-axis spin configuration. The Mn spin 
is experimentally known to be of the Heisenberg model 
with magnitude S' = 5/23 

The interaction between an itinerant spin and surroun- 
ding Mn spins in semiconducting MnTe is written as 

^i = -^J(r-^n)s-S„ (3) 

n 

where J(r — Rn) > is a ferromagnetic exchange inter- 
action between the itinerant spin s at r and the Mn spin 
Sn at the lattice site Rn • The sum on lattice spins 
is limited at some cut-off distance as will be discussed 
later. We suppose that J(r — Rn) is weak enough to be 
considered as a perturbation to the lattice Hamiltonian : 

J(r - Rn) = Ioexp[-a{r- Rn)] (4) 

where Iq and a are constants. We choose a = 1 for conve- 
nience. The choice of /q should be made so that the inter- 
action Hi yields an energy much smaller than the lattice 
energy due to H (see discussion on the choice of variables 
in Refs.^l^M). 
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Since in MnTe the carrier concentration is n = 4.3 x 
lO^^cm"^, low with respect to the concentration of its 
surrounding lattice spins o:^ lO^^cm"^, we do not take 
into account the interaction between itinerant spins. The 
contribution of that interaction to the energy of an itine- 
rant spin is negligible. 

We perform MC simulations on a thin film of dimen- 
sion L X L X L where L is the number of MnTe cells in 
y and z directions. Note that each cell contains two Mn 
atoms and two Te atoms (see Fig. [1]). Periodic boundary 
conditions are applied in all directions. The itinerant elec- 
trons move in the system under an electric field e applied 
along the x direction 

nE = -eel (5) 

where — e is the charge of electron and £ the displacement 
vector of the electron. Each electron spin s of magnitude 
1/2 interacts with neighboring lattice Mn spins within 
a sphere of radius Di according to Eq. (|3|). We take in 
the following the Ising model for the electron spin. In 
doing so, we neglect the quantum effects which are of 
course important at very low temperature but not in the 
transition region at room temperature where we focus 
our attention. 

The MC technique we use for the transport is a multi- 
step averaging which has shown to reduce efficiently sta- 
tistical fluctuations^^^^^. The reader is referred to our 
early works for a detailed description of our method as 
well as the effects of changing physical parameters such as 
Di in antiferromagnets. In a word, this technique consists 
in changing the lattice spin configuration as often as the 
temperature-dependent relaxation time allows. The re- 
laxation time of the lattice spin system is expressed as^^ 

" \1-T/Tn\'^ 

where A is a constant, u the correlation critical exponent, 
and z the dynamic exponent. From this expression, we 
see that as T tends to Tat, tl diverges. This phenomenon 
is known as the critical slowing-down. For the Heisenberg 
model, zu = 1.38 {u = 0.704 and z = 1.97).^^ We have 
previously shown that tl strongly affects the shape of 



III. RESULTS 

We have calculated the spin resistivity of the hexago- 
nal MnTe using the exchange integrals taken from Ref.^— . 
These values are given above. As said before, the values 
of the exchange interactions deduced from experimental 
data depend on the model Hamiltonian, in particular the 
spin model, as well as the approximations. Furthermore, 
in semiconductors, the carrier concentration is a function 
of T. In our model, there is however no interaction bet- 
ween itinerant spins. Therefore, the number of itinerant 



spins used in the simulation is important only for statis- 
tical average : the larger the number of itinerant spins 
the better the statistical average. The current obtained 
is proportional to the number of itinerant spins but there 
are no extra physical effects. 

The result of the spin resistivity is shown in Fig. [2j 
Using the value of p, we obtain the relaxation time of 
itinerant spin equal to r ~ 0.1 ps, and the mean free 
path equal to / c:^ 20A, at the critical temperature. 

P (Q.cm)^ 
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Figure 2: Spin resistivity p versus temperature T. Black 
circles are from Monte Carlo simulation, white circles are ex- 
perimental data taken from He et al.^^ . The parameters used 
in the simulation are Jl = — 21.5K, J2 = 2.55 K, J3 = —9 K, 
Io = 2K, D = 0.12 K, Di = a = 4.148A, ^ = 2 * 10^ V/m, 
L = 30. 

Several remarks are in order : 

i) With J3 slightly larger in magnitude than the value 
deduced from experiments, we find Tn = 310 K 

ii) p shows a pronounced peak in excellent agreement 
with experiments 

iii) Note that the shape of the peak depends on A. The 
value we used to obtain that agreement is A = 1 

iv) In the temperature regions below T < 140 K and 
above T/v the MC result is in excellent agreement with ex- 
periment, unlike in our previous work^^ using the Boltz- 
mann's equation 

v) In the region 140 K < T < Tat the MC result of p is 
slightly smaller than experimental data. In the search for 
an explanation, we see that the magnetization obtained 
by the MC simulation, though in good agreement with 
experiments at low T and yielding the precise value of 
Tat, is slightly smaller than the experimental one in the 
intermediate temperature region. This is shown in Fig. 
3. The magnetization deficit may be due to the fact that 
the magnetic anisotropy was too small D = 0.12 K ta- 
ken from Ref.^^, it is not strong enough to enhance the 
magnetization. 



IV. CONCLUSION 

We have shown in this paper the MC results of the 
spin resistivity p as a function of temperature in MnTe. 
We have taken into account the main interaction which 
governs the resistivity behavior, namely the interaction 
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between itinerant spins and the lattice Mn spins. Our re- 
sult is in agreement with experiments : it reproduces the 
correct Neel temperature as well as the shape of the peak 
at the phase transition. Note that the theory by Haas^^- 
predicts the absence of a peak in p in the temperature re- 
gion of the phase transition for antiferromagnetic MnTe. 
We finally emphasize that our excellent agreement was 
possible because most importantly we have correctly ta- 
ken into account the temperature dependence of the lat- 
tice spin relaxation time in the simulation^^i 



Figure 3: Magnetization M (in unit of /xs) versus tempe- 
rature T (in unit of K). Black circles are results of Monte 
Carlo simulation, white circles are experimental data taken 
from Efrem D'Sa et al.^°. The parameters are the same as 
those given in the caption of Fig. 2. 
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